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ABSTRACT 

We exploit a large, complete optical diameter and HI flux limited sample of spiral galaxies with 
types later than SOa to derive a robust measurement of the HI mass function (HIMF) for masses 
\og{Mni /Mq) > 7.4 which takes into account the effects of local large scale structure. The global 
HIMF derived for this optically-selected sample is well fit by a Schechter function with a = —1.24, 
log(M*/M ) = 9.99, 0» = 3.2 x 10" 3 Mpc -3 . These values match those derived from blind HI 
surveys to within the estimated uncertainties, yet our estimated HIMF is clearly lower than most 
other estimates at the lowest masses. We also investigate the variation in the derived HIMF among 
spiral subclasses, finding a clear distinction between the Schechter parameters found for types Sa-Sc 
and those Scd and later, in the sense that the HIMF of the latest types is rising at the low mass 
end, whereas that of the main spiral classes is flat or even declining. We also explore the possible 
environmental dependence of the HIMF by computing it separately in regimes of differing cosmic 
density. The HIMFs of higher density regions are found to have flatter low-mass ends and lower values 
of M* than those of lower density regions, although the statistical significance of the difference is low. 
Because the subsamples found in different density regimes exhibit virtually the same morphological 
fractions, the environmental dependence cannot be accounted for by morphological segregation, and 
must be a consequence of differences among galaxies of the same morphological type but found in 
different environments. If this dependence is caused by the well known deficiency of galaxies in clusters, 
then it would suggest that galaxies of small linear optical diameter are characterized by higher HI 
deficiency, an expectation consistent with gas removal mechanisms such as ram pressure stripping. 
Subject headings: galaxies: luminosity function, mass function — galaxies: distances and redshifts — 
radio lines: galaxies 



1. INTRODUCTION 

Because of its abundance in the interstellar medium 
and its emission strength at centimeter wavelengths, the 
neutral hydrogen content of galaxies serves as a good 
probe of galaxy environment and as an indicator of the 
potential for future star formation. For this reason, the 
probability distribution over HI mass, the HI mass func- 
tion (HIMF), can provide insight into the processes of 
galaxy formation and evolution. It is well known that 
a significant fraction of cluster galaxies are deficient in 
HI relative to field galaxies (see, e.g., Haynes, Giovanelli, 
& Chincarini 1984). Yet, to date, studies of the envi- 
ronmental dependence of the HIMF have been limited to 
comparisons of the HIMF for galaxies in or near a partic- 
ular cluster or void and those in the field (Briggs & Rao 
1993; Szomoru et al. 1996; Verheijen et al. 2001; Rosen- 
berg & Schneider 2002, hereafter RS; Davies et al. 2004; 
Gavazzi et al. 2004) . A more comprehensive study of the 
link between the environmental dependence of the HIMF 
and the deficiency of HI in clusters may shed light on the 
physical processes responsible for the gas depletion. 

While most early studies of the HIMF relied on 
optically-selected samples (e.g., Briggs & Rao 1993; 
Solanes, Giovanelli, & Haynes 1996, hereafter SGH), at- 
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tention has shifted more recently to results obtained from 
"blind" HI surveys (e.g., Zwaan et al. 1997; RS; Zwaan 
et al. 2003, hereafter Z03). Unlike optically-selected 
samples, Hl-selected samples have no bias against the 
low luminosity and low surface brightness galaxies which 
may be underrepresented in optical catalogs. However, 
as noted by SGH, HI blind surveys include their own se- 
lection effects. Notably, HI selected samples to date have 
tended to sample a much shallower volume than optical 
surveys, potentially leading to systematic errors in HI 
mass estimates due to distance uncertainties (Masters, 
Haynes, & Giovanelli 2004), and resulting in the selec- 
tion of a population of objects that is not representative 
of the population at larger distances. 

The largest HI blind survey to date is the HI Parkes 
All Sky Survey (HIPASS; Meyer et al. 2004, Zwaan et 
al. 2004). Z03 derived an HIMF from the HIPASS Bright 
Galaxy Catalog (BGC; Koribalski et al. 2004). Those 
authors also cross-correlated the BGC with known opti- 
cal galaxies, and identified previously-catalogued optical 
counterparts for 892 of the 1000 HI detections. As also 
found by RS, the HI blind samples include significant 
proportions of late-type low optical surface brightness 
systems that are not found in optical apparent magni- 
tude or diameter limited catalogs but are not otherwise 
unusual objects. The large fraction of BGC detections 
found to have previously-known optical counterparts was 
used to argue that, at least for the range of masses probed 
by the BGC, there is no evidence for a population of in- 
tergalactic HI clouds without stellar components (Korib- 
alski et al. 2004). 

Despite the focus on HI blind surveys, a large body of 
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HI data for optically selected galaxies already exists. If 
all (or nearly all) objects in the BGC mass range do in- 
deed have optical counterparts, then a sufficiently deep 
HI sample of optically selected targets may be able to 
sample a population of objects similar to that of the 
BGC. In particular, a sample that includes more sources 
and covers a greater volume than the BGC could allow 
for an investigation of the dependence of the HIMF on 
environment, as well as an independent check on the mor- 
phological dependence of the HIMF reported by Z03. 

The first two investigations of the HIMF based on HI 
observations of a large sample of optically selected galax- 
ies were Briggs & Rao (1993) and SGH. However, the 
former (based on 355 galaxies drawn from the catalog of 
Fisher & Tully 1981) is limited to a very shallow volume 
(cz < 1000 km/s). The latter (based on 532 galaxies 
drawn from the Arecibo General Catalog) was limited 
to morphological types Sa-Sc, and had significant uncer- 
tainty in the estimate of the low-mass end of the HIMF, 
owing to the fact that the sample did not include any 
objects with M HI < 10 8 5 M o . 

In this paper we present an HIMF derived from an 
optically selected sample of spirals and irregulars with 
both optical diameter and HI flux limits. The sample 
contains more than 2.7 times as many galaxies as any 
other used to derive an HIMF to date, and occupies a 
volume that greatly exceeds that of the BGC. We do 
not sample as deeply in HI mass however. We sample 
down to Mhi = 10 7A Mq , whereas Z03 reaches just be- 
low Mhi = 10 7 -°M Q . Because of the large volume sam- 
pled, we also investigate the dependence of the HIMF on 
morphology and environment, and compare our results 
to those of other authors. 

2. THE DATA 

2.1. Sample Selection 

Like SGH, our sample is drawn from the all-sky pri- 
vate database maintained by M.P.H. and R.G. known as 
the Arecibo General Catalog (AGC). M.P.H. and R.G. 
also maintain a digital archive of ~9100 HI line spectra 
of AGC galaxies, first described in Haynes et al.{ 1999). 
The majority of the galaxies in the archive are in the 
'Arecibo sky' (-2° < deck < +38°). HI observa- 
tional data for galaxies in that region was presented in 
Giovanelli & Haynes (1985a), Giovanclli et al. (1986), 
Giovanelli & Haynes (1989), Wegner, Haynes, & Gio- 
vanclli (1993), Giovanelli & Haynes (1993), Giovanelli et 
al. (1997), Haynes et al. (1997), and Haynes et al. (1999) 
and also includes new observations of 120 galaxies made 
with the Arecibo 305 m telescope of the National As- 
tronomy and Ionosphere Center. All of these data have 
been newly reprocessed, using the HI flux and width cor- 
rections described in Haynes et al. (1999). This homoge- 
neous, archival dataset is being prepared for publication 
(Springob, Haynes, & Giovanelli 2004, in preparation) 
and access to the digital spectral archive will be pro- 
vided through the U.S. National Virtual Observatory. 
As in Haynes et al. (1999), we augment the sample of 
galaxies in the archive with HI fluxes from the literature 
as compiled in the AGC. The total number of HI flux 
measurements available in this sample is 13,773. 

Both methods used to derive the HIMF (described in 
Section 3) require an accurate characterization of the se- 
lection criteria of the galaxy sample. We thus need to 



select a subset of the galaxies from the archive which 
is complete for a particular set of criteria. The AGC 
contains all objects in the Uppsala General Catalog of 
Galaxies (hereafter UGC, Nilson 1973), which claims 
to be complete down to an apparent diameter of 1.0' 
for decl. > —2.5°. Thus, almost all AGC galaxies in 
the Arecibo declination range with apparent diameter 
greater than 1.0' are UGC objects. There are a few ex- 
ceptions to this, as the UGC is not quite complete close 
to this diameter limit. In some instances, double galax- 
ies listed as single in the UGC are listed singly in the 
AGC, and some additional objects, notably those drawn 
from the Catalog of Galaxies and Clusters of Galaxies 
(Zwicky et al. 1960-1968) that are large enough (a > 1.0') 
are included. The diameter completeness of the UGC is 
discussed in detail in Hudson & Lynden-Bell (1991). 

We therefore select a sample of 2771 galaxies (hereafter 
referred to as the "complete sample" , as distinct from the 
subsamplcs discussed in Sections 4.1 and 4.2), containing 
all objects in the HI digital archive or the AGC with 
a > 1.0 (where a is the optical apparent diameter in 
arcminutcs), —2° < deck < +38°, galactic latitude \b\ > 
15°, morphological type Sa-Irr, and for which log(FHi) > 
0.6 (where Fhi is the corrected HI flux density integrated 
over the profile in units of Jy km/s). Ellipticals, SOs, and 
SOas are excluded due to the poor completeness of the 
archive for these types (see Section 2.3). The choice of 
HI flux limit is explained in Section 2.3. 

2.2. Galaxy properties 

2.2.1. Optical properties 

Morphological types and optical diameters were de- 
rived using the same methods described in Section 2.2 
of SGH. As in SGH, we define the linear optical diame- 
ter, D 0} as the physical size of major diameter a as seen 
at a radial distance r. This is given by the relation 

(ZVkpc) = 0.291(r)a (1) 

where r is expressed in units of Mpc. A Hubble constant 
of H n = 70 km/s Mpc -1 is assumed here and throughout 
this paper. 

2.2.2. Radial distance 

Heliocentric radial velocities are derived from the HI 
line spectrum. The distribution of the systemic velocities 
for the complete sample is shown in the upper left panel 
of Figure 1 . It reflects strongly the principal features in 
the local large scale structure in this area of sky, most 
notably the Local, Pisces-Perseus and Coma-A1367 su- 
perclusters. As Masters, Haynes, & Giovanelli (2004) 
point out, however, local departures from Hubble flow 
are likely to be significant enough that we cannot simply 
use the redshift distance for each galaxy. We thus con- 
vert the velocities into distances using the multi-attractor 
flow model of Tonry et al. (2000), modified to assume our 
value of Ho — 70 km/s Mpc -1 . The model includes in- 
fall onto Virgo and the Great Attractor, but reduces to 
v/H at large distances. 

For galaxies in clusters (or large groups), however, the 
peculiar velocity induced by the cluster or group is likely 
to make the Tonry model applied to the galaxy redshift 
an insufficient estimate of the distance. Thus, where pos- 
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Fig. 1. — Distribution of heliocentric radial velocities for all 
galaxies in both the complete sample and the n < 1.5, 1.5 < n < 
3.0, and n > 3.0 subsamples, in bins of width 200 km/s. The 
peak velocity of the complete sample at ~ 4800 km/s is roughly 
a factor of three greater than that of the BGC (see Koribalski et 
al. 2004), the largest HI selected sample for which a HIMF has 
been derived. The mean velocity is similar to that of the Arccibo 
Dual-Beam Survey used by RS (see Rosenberg & Schneider 2000 
Figure 12), but our sample provides much better statistics. Note 
how the denser subsamples are "clumpier" than the lowest density 
regions. 



sible, we assign group or cluster membership to galax- 
ies and use the group or cluster redshift instead. 68 of 
the galaxies in our sample have been identified as clus- 
ter galaxies by either Giovanelli et al. (1997) or Dale et 
al. (1999). For each of these objects, we use the cluster 
redshift. An additional 777 galaxies are either identified 
in the "P2" Nearby Optical Galaxies (NOG) group cata- 
log (Giuricin et al. 2000), or were found to be NOG group 
members via a "friends-of-friends" algorithm similar to 
that used by the P2 catalog. For these objects, we use 
the NOG group redshift. 

2.2.3. HI Mass 
The neutral hydrogen mass of each source is given by 



MhiIMq = 2.36 x 10 VF ffl (2) 

HI fluxes are corrected for beam attenuation, HI self- 
absorption, and pointing offsets as described in Springob, 
Haynes, & Giovanelli (2004, in preparation). 

2.3. Completeness 

The archive contains HI data for 74% of Sa-Sb galaxies, 
88% of Sbc-Sc galaxies, 81% of Scd-Sd galaxies, and 97% 
of Sm-Irr AGC galaxies with diameter a > 1.0 ' in the 
sky region -2° < dccl. < +38°, \b\ > 15°. In contrast, 
it contains a mere 14% of E-SOa galaxies fulfilling these 
criteria. These fractions include nondetections and obvi- 
ous cases of HI confusion. However, as shown in Figure 
2, virtually all of the galaxies without HI data lie close 
to the diameter limit. Imposing a suitable HI flux limit 
should improve the completeness considerably. 

We can crudely estimate the HI flux completeness of 
the sample by relying on the linear relationship between 
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Fig. 2. — Diameter completeness for different morphological sub- 
samples. Thin line shows distribution of apparent diameters for all 
galaxies with a > 1.0 in the AGC in the sky region —2° < decl. 
< +38° and |6| > 15°. Thick line shows the number of those galax- 
ies for which there arc HI parameters in the HI archive. Bins have 
width 0.1 dex. The distributions do not show any noticeable devia- 
tions from the overall diameter distribution for the UGC presented 
in Hudson & Lynden-Bell 1991. 



log(_D ) and log(M#j) given by Haynes & Giovanelli 
(1984). This relationship is plotted for our sample in Fig- 
ure 3. A linear least-squares fit gives log(Z?o/kpc) = 0.46 
\og(Mni /M C3 ) — 3.06, with an average scatter of 0.12. 
Since the scatter is nearly Gaussian, we can estimate, 
for any galaxy of HI mass Mhi, the fractional number 
of galaxies with that mass, at that distance, which are 
not included in the sample because they fall below the 
diameter limit. From this we can infer a "corrected flux 
distribution" (Figure 4), the distribution of the HI fluxes 
that we would have if there were no diameter limit. For 
an HI flux-limited sample, the logarithm of the number 
of galaxies per unit log(Fui) should be proportional to 
-1.5 log(Fff/). As Figure 4 shows, if one normalizes this 
relation to include the same number of galaxies as the 
actual data, the completeness drops to less than 80% per 
bin below log(Fni) = 0.6, so we adopt this as our HI 
flux limit. In contrast, SGH actually used a lower flux 
limit of log(Fjji) = 0.4, but did not reach nearly as far 
down in HI mass due to their reduced sky coverage and 
restriction to types Sa-Sc. Z03 used a peak flux limit, 
rather than an integrated flux limit as we have done, of 
116 mJy As seen in Figure 9 of Koribalski et al. 2004, 
this means that all but a handful of BGC objects have 
integrated fluxes of at least two or three times greater 
than our limit of 10 6 Jy km/s. 

The AGC includes 107 galaxies with a > 1.0, -2° < 
decl. < +38°, galactic latitude |6| > 15°, and type Sa-Irr 
for which HI observations yielded nondetections. How- 
ever, only 7 of these galaxies have an rms noise per chan- 
nel greater than 5 mJy. If we were to assume that the 
emission profiles for these galaxies are rectangular with 
amplitude 1.5 times the rms noise per channel, then a 
galaxy with rms noise equal to 5 mJy would require a 
velocity width of at least 530 km/s to be included in a 
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Fig. 3. — Logarithm of linear optical diameter, Do, vs. logarithm 
of HI mass, Mhi- Crosses give the values for all galaxies in our 
complete sample. Open circles are the values for the BGC for 
which optical diameters can be found in the AGC. Solid line shows 
best-fit linear relationship between Do and Mhi for our sample 
only, log(A)/ kpc) = 0.46 log(M HI /M Q ) - 3.06. 



Fig. 4. — Logarithmic "corrected flux distribution". Filled cir- 
cles show the number of galaxies per 0.1 dex bin in our sample, 
corrected to include galaxies that are missed due to the diameter 
limit. A flux limited sample should give <i(log N ga [ axies )/d(\og 
Fhi) —1.5 log Fhi + constant. This solid line shows this re- 
lation, normalized to include an equal number of objects as the 
actual data for log Fhi > 0.6. 



flux limited sample with flux limit 10°' 6 Jy km/s. It thus 
seems exceedingly unlikely that the nondetections would 
have much impact on our result, so we exclude them from 
the derivation of the HIMF explained in Section 3. 

There are also 55 galaxies that would otherwise meet 
our selection criteria, but are excluded because they are 
clearly confused, and it is impossible to determine how 
much of the HI flux belongs to the target galaxy. Both 
the redshift and optical diameter distributions of these 
galaxies seem to match that of the complete sample, so 



we do not believe the exclusion of these galaxies has a 
significant effect on the HIMF shape. 

3. THE HI MASS FUNCTION 
3.1. Methods 



3.1.1. 



Method 



The T,(l/V ma x) method, first developed by Schmidt 
(1968), has frequently been applied to the determination 
of HI mass functions (e.g., Zwaan et al. 1997; RS). Each 
galaxy is assigned a weight given by 1/V max (xi, x 2 , ...), 
where V max (xi, x 2 , • ■•) is the effective search volume for 
a galaxy with a particular set of physical properties 
Xi,X2---- relevant to the observational selection criteria. 
The HIMF, 4>(Mhi), is defined as the space density of 
objects with mass Mhi- And so the value of ()){Mhi) for 
a particular mass bin is given by the sum of the weights 
of all galaxies in that bin. 

In our case, the selection function depends on optical 
diameter and HI flux. Following equations 1 and 2, we 
see that 



V max (D ,M HI ) = ^h*Ly V un 
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M Q 2.36 x 10 5 F Hmm 
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where an 



, llm and FHiiim are the optical diameter and HI 
flux limits (1' and 10 6 Jy km/s respectively), and f 8 k y 
is the fraction of the sky covered by the survey (equal to 
0.257 for -2° < decl. < 38° and \b\ > 15°). 

The primary concern about the T,(l/V max ) method 
is the implicit assumption that sources are distributed 
homogeneously within the survey volume. The survey 
volume for our sample contains rather significant den- 
sity inhomogeneities, most notably, the Local Superclus- 
ter and the Pisces-Perseus Supercluster, so noticeable 
in Figure 1. To offset this effect, modified versions of 
T,(l/V max ) which account for the effects of large scale 
structure have been used by RS and Masters, Haynes, 
& Giovanelli (2004). In both of these cases, each galaxy 
was weighted by l/V max - e ff ective = l/[n(V max )V max )} 
rather than 1/V max , where n(V max ) is the average den- 
sity within the volume V max normalized by the average 
density of the universe. 

However, while RS derived n{V max ) from a magnitude- 
limited sample of galaxies drawn from the NASA/IPAC 
Extragalactic Database, Masters, Haynes, & Giovanelli 
(2004) used the matter density reconstruction derived 
from the IRAS Point Source Catalog Redshift Survey 
(PSCz). The latter reduces the likelihood that the esti- 
mate of n(V max ) will be influenced by distortions in red- 
shift space due to galaxy peculiar velocities, since such 
peculiar velocities are accounted for by the density recon- 
struction. We have therefore chosen to use the Masters, 
Haynes, & Giovanelli (2004) approach. 

The PSCz density reconstruction (presented in Bran- 
chini et al. 1999) consists of a cartesian grid of mass den- 
sities (with grid spacings of 1.34 Mpc out to a limit of 
85.7 Mpc) derived from the distribution of PSCz sources 
in redshift space, assuming (3 = £l^ as8 /b = 0.5. At radial 
distance steps of 1.34 Mpc, we average all PSCz density 
gridpoints interior to the radial distance at the step that 
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occupy -2° < decl. < 38° and |b| > 15°. The average 
density for that volume is taken to be n(V m ax)- Beyond 
85.7 Mpc, we assume n(V max ) = 1. 

3.1.2. 2D S WML Method 

An alternative method for deriving the HIMF 
is the two-dimensional stepwise maximum likelihood 
(2DSWML) method, developed by Loveday (2000) for 
constructing an optical luminosity function, and first 
used for the HIMF by Z03. Z03 used this method rather 
than T,(l/V max ) because of the latter's insensitivity to 
the effects of large scale structure. We believe this con- 
cern has been abrogated by our use of a modified version 
of T,(l/V max ) that accounts for large scale structure as 
described in the previous subsection. 2DSWML has its 
own drawbacks: Unlike T,(l/V max ), it is not automati- 
cally normalized and it also assumes that the shape of 
the mass function is the same everywhere. The latter 
effect can skew the HIMF at the high-mass end to match 
this assumption. 

Nevertheless, we have implemented the 2DSWML 
method as a check on our result. Rather than binning 
the sample by HI mass and HI velocity width, as Z03 
have done, we have binned by HI mass and linear optical 
diameter. Bins are 0.25 dex wide in logarithmic mass 
(just as with the T,(l/V max ) method) and 0.1 dex wide 
in logarithmic diameter. A detailed description of the 
method is provided by Z03. 

4. RESULTS 

Figure 5 shows the HIMF, as determined by both the 
T,(l/V max ) and 2D S WML methods, for the complete 
sample. The normalization for the 2DSWML points 
has been chosen to minimize scatter with the Yj(l/V max ) 
points. Agreement between the two methods appears to 
be quite good, with the 2DSWML points falling within 
lcr of the T,(l/V max ) for 12 of the 15 bins (and within 
0.2cr for 7 of those bins). We note, however, that the 
high-mass end of the 2D S WML HIMF is influenced by 
values in lower mass bins (as mentioned in Section 3.1.2). 
The two bins that show the greatest deviation from the 
^0-/V max ) results are in fact the two highest mass bins. 
Since the earliest types tend to have very little HI (see, 
e.g., Roberts & Haynes 1994), we expect that the exclu- 
sion of types earlier than Sa makes little difference in the 
resulting HIMF. 

We have fit the T,(l/V max ) HIMF to a Schechter func- 
tion (Schechter 1976) defined by 

d \og{M HI ) 

= ln(10)<MMWM«r +1 cxp(-MWM») (4) 

Our best fit Schechter parameters are a = —1.24, 
log(M*/M ) = 9.99, c/>» = 3.2 x 10" 3 Mpc" 3 , with a 
xt of 0.99. Our Yj(l/V max ) errorbars are based on Pois- 
son counting statistics, and ignore any errors introduced 
by binning. We have experimented with varying the bin 
sizes across the range of Mhi (so that the bins are more 
finely spaced for higher values of Mhi)- This does yield 
a modest decrease in ^i/ 2 , but no significant change in 
the values of the Schechter parameters, so we stick with 
the fixed bin widths for the sake of simplicity. 
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Fig. 5. — Top, HIMF for the complete sample, using the 
^(1/Vmax) method (filled circles) and the 2DSWML method (open 
squares). T,(l/V m ax) errorbars indicate lcr uncertainties based 
on Poisson counting statistics. The thick solid line is the best- 
fit Schechter function to the T,(l/Vmax) points: a = —1.24, 
log(M*/M ) = 9.99, <t>* = 3.2 X 1CT 3 Mpc" 3 . Other lines show 
the Schechter parameters from five HI blind surveys. (Zwaan et 
al. 1997: thick dotted; RS: thin dashed; Henning et al. 2000: thick 
dashed; Kilborn 2000: thin dotted; Z03: thin solid); bottom, distri- 
bution of HI masses in the full sample per 0.25 dex bin. Comparing 
this distribution to Rosenberg & Schneider (2000) Figure 16 and 
Z03 Figure 4, one can sec that, while we have many more objects 
across most of the mass range, we do not go any lower in mass than 
those surveys, and in fact have fewer objects at the lowest masses 
than Z03. 



Best fit values for a and M#, along with error con- 
tours (based on x 2 goodness of fit), are shown in Figure 
6. For comparison we plot the results for the five HI se- 
lected samples listed in Tabic 3 of Z03 (i.e., the results of 
Zwaan et al. 1997, Kilborn 2000, Henning et al. 2000, and 
RS, as well as Z03 itself). The agreement between our 
Schechter parameters and those of HI selected samples 
is no worse than the agreement of the latter's Schechter 
parameters with one another. As Masters, Haynes, & 
Giovanelli (2004) point out, however, the slope of the Z03 
HIMF is actually steeper when one takes into account the 
effects of peculiar velocities. The authors estimate that 
a flow model-corrected Z03 HIMF would actually have 
a low mass slope of -1.4, which would put our result on 
the "shallow" end of HIMF results along with Zwaan et 
al. (1997). In fact, even if we accept the Z03 HIMF as 
published, one must account for the authors' high value 
for the normalization constant 0* in making any com- 
parison of the HIMFs at the low mass end. In Figure 5, 
we have overplotted the HIMFs of all five HI blind sam- 
ples using the Schechter parameters as published. At 
masses of ~ 10 8 M Q , both our HIMF and that of Zwaan 
et al. (1997) are lower than each of the others by a fac- 
tor of a few. Another way to compare the results is by 
comparing the total HI mass density, derived by inte- 
grating the Schechter function over all HI masses. Using 
Phi = <t>*T(2+a)M*, we get p m = 3.8 x 1O 7 M /Mpc" 3 . 
Converting this to Q H i gives Sl H i = 2.7 x 10~ 4 , 39% 
lower than the value reported by Z03. 

There are several possible explanations for this dis- 
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TABLE 1 

SCHECHTER PARAMETERS 




log(M./M e ) 

Fig. 6.— The \ 2 contours of the Schechter fit to the HIMF 
for the complete sample. Contours represent 1, 2, and 3<r. Also 
shown are Schechter parameters from five HI blind surveys: the 
Arecibo HI Strip Survey (AHISS, Zwaan et al. 1997), the Arecibo 
Dual-Beam Survey (ADBS, RS), the HIPASS HIZSS (Henning et 
al. 2000), the HIPASS SCC (Kilborn 2000), and the HIPASS Bright 
Galaxy Catalog (BGC, Z03). ADBS crrorbars are estimated from 
the contours in RS Figure 3. AHISS crrorbars are from Zwaan 
(2004, private communication). 



crepancy. It may be that galaxies with small physical 
diameters, not included in our complete sample because 
of our diameter limit, are large contributors to the low 
mass end of the HIMF. In Figure 3, we have overplot- 
ted log(Z)o) vs. log(Mjyj) for all BGC galaxies for which 
there are optical diameters listed in the AGC. As with 
our own data, the values of Dq and Mhi are calculated 
using the Tonry model-corrected distance. It does ap- 
pear that across most of the Mhi distribution, the BGC 
includes a population of galaxies with smaller diameters 
than those included in our complete sample. However, 
this is least apparent for log(M^//M Q ) < 8.5, where the 
divergence between the HIMFs is greatest. 

We also note that Briggs & Rao (1993), the only other 
work to estimate the HIMF based on a large optically se- 
lected sample including the range of morphological types 
we include here, finds a low mass slope of a ~ —1.25, 
nearly identical to our value. This suggests that the dis- 
crepancy may be due to some hidden systematic bias 
in the selection of optical samples. HI selected samples 
may be detecting faint low surface brightness galaxies 
that are missed by optical surveys. This could include 
galaxies that meet our diameter limit, but arc nonethe- 
less not included in our sample because they are too faint 
to be optically detected. Yet we find that, of the 758 
BGC galaxies with \b\ > 15°, only 39 do not have optical 
counterparts with diameter estimates listed in the AGC. 
And only 3 of those have log(M HI /M e ) < 8.5. 

Another possibile explanation for the disagreement be- 
tween our estimate of the low mass end of the HIMF 
and those of most HI selected samples is that the dis- 
crepancy is a consequence of the fact that the lowest 
mass galaxies in each of these samples are confined to a 
very local volume. In both the BGC and our complete 
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sample, all galaxies with log(M#j/M Q ) < 8.5 are found 
within Vheiio < 1200 km/s. There could be systematic 
distance errors due to defects in the flow model, or there 
could be systematic environmental effects within that lo- 
cal volume. We note that the only HIMF derived from an 
HI blind survey which matches our result at low masses 
is Zwaan et al. (1997), which, like our sample and un- 
like HIPASS, includes nearby galaxies in the direction of 
Virgo. 

It should be noted that the Zwaan et al. (1997) sam- 
ple, like the majority of ours, was observed with the 
305m Arecibo telescope, which has a much narrower 
beam (3.3') than the Parkes 64m radio telescope used for 
HIPASS (15.5'). Additionally, with few exceptions, in- 
teracting and highly disturbed systems are not included 
in our sample. We thus do not believe that confusion 
has any meaningful impact on the HIMF shape we have 
derived for our sample, but it may have some impact on 
the HIMFs derived from HIPASS. It is not clear whether 
confusion would systematically effect low or high mass 
objects however. 

Given both the distance uncertainties and our poor 
understanding of the environmental dependence of the 
HIMF (which we discuss in Sections 4.2 and 5), we can- 
not resolve the question of why our HIMF diverges from 
that of most HI blind surveys at the low mass end. These 
effects can only be disentangled with the aid of a deeper 
sample that includes low mass galaxies at larger dis- 
tances. 

4.1. Morphological Dependence 

We have computed the HIMF for four subsets of the 
complete sample, divided by morphological type: Sa-Sb, 
Sbc-Sc, Scd-Sd, and Sm-Irr. The results are plotted in 
Figure 7, and best-fit Schechter parameters are given in 
Table 1. Figure 8 shows the 1 a contour plots, along 
with the corresponding results from Z03 and SGH. (Our 
morphological types were derived in the same manner as 
in SGH, i.e., primarily using UGC types. Z03, on the 
other hand, uses Lyon/Meudon Extragalactic Database 
types.) Our results are broadly consistent with those au- 
thors for most morphological types. Types Sa-Sc occupy 
a parameter space that is clearly distinct from types Scd 
and later: The former grouping has higher values of both 
M* and a than the latter. 

Z03 makes no estimate of the uncertainty on their best 
fit Schechter parameters for their morphological subsam- 
ples. However, it appears that if their uncertainties are 



7 




log(M HI /M s ) 



Fig. 7. — HIMFs for the each of the morphological subsamplc 
Solid lines show best-fit Schcchtcr functions. 
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Fig. 8.— The \ 2 contours of the Schechter fit to the HIMF 
for morphological subsamples Sa-Sb, Sbc-Sc, Scd-Sd, and Sm-Irr. 
Contours represent 1 a uncertainties. Corresponding parameters 
from Z03 and SGH arc also shown. 



similar in magnitude to ours, then the 1 a uncertainties 
would overlap or nearly overlap for all types with the 
possible exception of Scd-Sd. The discrepancy for that 
case could be due to differences in the morphological type 
assignment criteria. 

The higher values of both M* and a for Sa-Sc galaxies 
than for later types is not unexpected, given the decrease 
in the average HI mass for later types relative to Sa-Sc 
shown in Table 1 of Roberts & Haynes (1994). 

4.2. Environmental Dependence 

We have also computed the HIMF for three subsets 
of the complete sample, divided by local matter density 
as determined by PSCz (n < 1.5, 1.5 < n < 3.0, and 
n > 3.0). For this task, the version of the PSCz density 
reconstruction referred to in Section 3.1.1 is insufficient, 



due to its 85.7 Mpc limit. We therefore use an alternate 
version of the density reconstruction, with grid spacings 
of 2.68 Mpc out to a limit of 171.4 Mpc. 

In order to compute the HIMF for each of these subsets 
using the T,(l/V max ) method, we must account for the 
fact that, not only does the space density of objects vary 
with distance due to large scale structure, but the frac- 
tional volume of space occupied by regions of a particular 
density does as well (see Figure 1). We thus compute the 
HIMF for each of our three subsamples using the effective 

V m a.r. 



^max — effective — TYlifl 

fint {Vmax ) 
^int (Vmax) 



(171.4 Mpc) 2 



(5) 



where V max is given by Equation 3, fi n t(V max ) is the 
fraction of the volume V max occupied by regions of the 
density range for each given subsample, and ni nt (V max ) 
is the average density within those regions. n(V max ) now 
represents the average density within regions of the den- 
sity range for each given subsample, not within the entire 
volume V max . Since PSCz only extends to 171.4 Mpc, all 
galaxies with distances greater than 171.4 Mpc are dis- 
carded, and V max - e f fective is capped at that distance. 

The resulting HIMFs are shown in Figure 9, with 
corresponding \ 2 contour plots in Figure 10. Best-fit 
Schcchtcr parameters can be found in Table 1. We find 
that the lowest density sample has a steeper low-mass 
slope and a higher value of than the two higher den- 
sity samples. This trend is consistent with the results of 
Briggs & Rao (1993), RS, and Davies et al. (2004), all of 
whom found evidence for a flatter HIMF in Virgo than 
in the field, and Verhcijen et al. (2001), who found a flat- 
ter HIMF in Ursa Major than in the field. However, as 
seen by the overlapping 1 a error contours, the uncertain- 
ties are large enough that the data are consistent with 
the Schechter parameters being identical for all three re- 
gions. The data are insufficient to show conclusively that 
the HIMF varies with environment. 

We also note that, because the density field is derived 
from a flux limited catalog (and smoothed with a Gaus- 
sian filter of constant radius 4.57 Mpc) in which the av- 
erage galaxy-galaxy separation increases with distance, 
the shot noise of the density field increases with distance. 
Thus, for example, the set of all gridpoints with observed 
n > 3.0 at large distances will contain more contami- 
nation from less dense regions than the corresponding 
set of gridpoints with observed n > 3.0 at nearby dis- 
tances. This effect may have some influence on our re- 
sults, though it appears not to create any deviation of the 
HIMF shape from a Schechter function. As seen in Table 
1, xt is actually lower for all three of the density-selected 
subsets than it is for the complete sample. 

Finally, we check whether the apparent change in the 
slope of the low mass end of the HIMF is a consequence 
of the correlation between a and M». If M* varies with 
environment, then it may result in a change in a that 
does not represent a true change in the low mass slope. 
We have computed the least-squares best fit line for all 
points with \og{M HI /M & ) < 9.5 in the log((f>)-log(M HI ) 
plots of Figure 9 for each of the three density regimes. 
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Fig. 9.— HIMFs for the n < 1.5 (open circles), 1.5 < n < 3.0 
(open squares), and n > 3.0 (open triangles) samples. Solid, long- 
dashed, and short-dashed lines show best-fit Schechter functions 
for each of the three samples respectively. 
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Fig. 10.— The \ 2 contours of the Schechter fit to the HIMF 
for the n < 1.5, 1.5 < n < 3.0, and n > 3.0 samples. Contours 
represent 1 a uncertainties. 



(Cowie & Songaila 1977), and turbulent viscous stripping 
(Nulsen 1982). Valluri & Jog (1991) showed that all of 
these mechanisms should be more pronounced for galax- 
ies with small optical sizes. Their data, however, showed 
that galaxies with smaller optical sizes tend to be less 
deficient in HI than galaxies with larger optical sizes. 
The authors used this result to rule out any mechanism 
for gas removal that involves galaxy-ICM interactions. 
The only mechanism investigated by Valluri & Jog (1991) 
that the authors concluded might be consistent with the 
diameter-dependent deficiency shown by their data was 
galaxy-galaxy tidal encounters in subclumps within the 
cluster. This interpretation was criticized by Solanes et 
al. (2001), who noted that the deficiency parameter used 
by Valluri & Jog (1991) had a residual dependence on 
the optical diameter. The data presented by Solanes et 
al. (2001) suggested that in fact, there is no apparent 
dependence of HI deficiency on optical diameter. 

If HI deficiency in clusters does depend on optical size, 
then the dependence should manifest itself in the cluster 
HIMF. //the population of galaxies in clusters was oth- 
erwise identical to the population of galaxies in the field, 
then a constant deficiency for all cluster galaxies would 
mean a cluster HIMF with the same shape as the field 
HIMF, but shifted to lower HI masses. 

Of course, not all cluster galaxies have the same HI 
deficiency, so we require a somewhat more sophisticated 
model. If we define 4> c i(Mhi) as the cluster HIMF that 
would exist if none of the galaxies were deficient, then we 
can express <fi' cl (MHi), the cluster HIMF that includes HI 
deficient galaxies, as 

<P'ci(Mhi) = MM hi )4> c1 (Mhi) + fiMMHi + Aj) 

(6) 

where fo(Mui) is the fractional number of galaxies with 
HI mass Mhi that are not HI deficient, fi is the fractional 
number of galaxies with expected HI mass Mhi + A, that 
are deficient by a factor Aj, and the summation is over 
all expected HI masses greater than Mhi- We define the 
deficiency DEF as in Haynes & Giovanelli (1984), the log- 
arithm of the ratio of the expected HI mass for a particu- 
lar linear optical diameter and morphological type to the 
observed HI mass for that galaxy. At the low mass end 
of the HIMF, where we can approximate log(<f> c i(MHi)) 
~ (a + l)log(M hi) +constant, (6) becomes 



We find slopes of -0.51, -0.32, and -0.31 for the n < 
1.5, 1.5 < n < 3.0, and n > 3.0 regions respectively, 
confirming that the steeper slope in the n < 1.5 regions 
is not an artifact of the fitting parameters. 

5. DISCUSSION 

Studies of the HI content of galaxies in clusters have es- 
tablished that a significant fraction of cluster galaxies are 
deficient in HI gas relative to the field (e.g., Giovanelli & 
Haynes 1985b). The most commonly suggested mecha- 
nisms invoked to explain the removal of gas from cluster 
galaxies involve one form or another of interaction be- 
tween the HI gas in galaxies and the intercluster medium 
(ICM). These mechanisms include ram-pressure strip- 
ping (Gunn & Gott 1972), thermal evaporative stripping 



^ci(Mhi) = (fo(MHi)+^2f^ (a+1)Al ) 

^ i ' 

^-^Q(a+l)log(MHi)+constant 

If /o and fi are completely independent of linear op- 
tical diameter (that is, independent of the expected HI 
mass), then deficiency will not affect the low mass slope. 
If, on the other hand, either the fractional number of 
HI deficient galaxies or the amount of HI deficiency per 
galaxy increases (decreases) with diameter, then the low 
mass slope will get steeper (shallower). 

None of this analysis necessarily applies to a realistic 
comparison of the cluster and field HIMF. Morphological 
segregation in high density regions means that the field 
HIMF cannot be represented by 4> c i{Mhi)- As shown in 
Section 4.1, the HIMF for early type spirals has a much 
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flatter low mass slope than the HIMF for later types. So 
a flatter HIMF in high density regions could be a natural 
consequence of early type spirals being more abundant 
than later types in these regions. 

Among the three levels of PSCz density that we study 
here, however, there is essentially no morphological seg- 
regation. 34% of the n < 1.5 galaxies are Sa-Sb, 45% 
are Sbc-Sc, 6% are Scd-Sd, and 15% are Sm-Irr. For 
the 1.5 < n < 3.0 and n > 3.0, the breakdowns are 
36%/46%/7%/12% and 37%/44%/7%/12% respectively. 
We have constructed synthetic HIMFs by adding linear 
combinations of the HIMFs for each morphological type, 
described in Section 4.1, in the ratios for which they ap- 
pear in each of these three density regions. We have 
then computed best-fit Schechter functions for each of 
these HIMFs using the same errorbars on each bin shown 
in Figure 9. The resulting low mass slopes are a = - 
1.27, -1.21, and -1.25 for the n < 1.5, 1.5 < n < 3.0, 
and n > 3.0 regions respectively, a much smaller devi- 
ation than is actually observed. So we conclude that if 
the marginally flatter HIMF in the two higher density 
regimes is a result of increased HI deficiency relative to 
the low density regime, then we believe this would be an 
indication of increasing deficiency towards smaller opti- 
cal diameters at the low mass end of the HIMF. 

As a quantitative measure of deficiency, we use the ex- 
pected D — Mfjj relations from SGH for types Sa-Sc 
and the relations from Haynes & Giovanelli (1984) for 
types Scd-Irr. The resulting distribution of deficiency as 
a function of optical diameter is shown for all three den- 
sity subsamples in Figure 11. The average deficiency for 
each subsample is slightly negative (~ —0.2 for all three 
subsamples). This is not surprising, given that HI defi- 
cient galaxies are more likely than non-deficient galaxies 
to fall below our HI flux limit and not make it into the 
sample. We can see from these plots that the fractional 
number of HI deficient galaxies (DEF> 0.3, as defined 
by Solanes et al. 2001) in the n < 1.5 and 1.5 < n < 3.0 
samples is relatively small (2% in both cases). In the 
n > 3.0 sample, the number is somewhat larger (8%), 
but the only apparent trend of HI deficiency is a some- 
what increased number of deficient galaxies towards in- 
termediate diameters. Of course, our selection criteria 
discriminate against HI deficient galaxies with small di- 
ameters, so there could very well be a large population 
of deficient galaxies with small diameters that is being 
missed by our sample. 

To answer this question, we examine the sample of all 
AGC galaxies of types Sa-Sc with optical apparent di- 
ameter greater than 1 arcminute that can be found in 
the sky region covered by our complete sample. For each 
of the three density regions, we compute the fraction of 
these galaxies with HI fluxes greater than our HI flux 
limit. In Figure 12, we show the ratio of these "HI de- 
tection fractions" between the different density regions, 
divided into bins of physical linear diameter. While the 
uncertainty in the ratio for each bin is significant, there is 
a clear trend in that the fraction of Hi-detected galaxies 
drops towards smaller diameters more sharply in both 
the intermediate and high density regions than in the 
lowest density regions. Since these plots are restricted 
to Sa-Sc and, as shown in Section 4.1, the Sa-Sb and 
Sbc-Sc subsamples have nearly identical HIMFs, it is un- 
likely that this is due to any morphological segregation 
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Fig. 11. — Deficiency vs. linear optical diameter for the n < 1.5, 
1.5 < n < 3.0, and n > 3.0 samples. 



Fig. 12. — Ratio of the fraction of Sa-Sc galaxies larger than 1 
arcminute in optical angular diameter that meet our HI flux limit in 
higher density regions to that same fraction in the n < 1.5 regions, 
plotted against physical optical diameter. The filled circles are for 
the n > 3.0 sample, and open circles are for the 1.5 < n < 3.0. 
While the errorbars are large, there is a clear trend towards lower 
HI dctcctability in smaller galaxies for the high and intermediate 
density regions than there is for the lowest density regions. 



effects. In fact, this may explain why we observe so little 
morphological segregation among our three environmen- 
tal subsamples. As noted in Solanes et al. (2001), HI 
deficiency is preferentially found in galaxies with earlier 
morphological types. While they may be more numerous 
in high density regions, they are also harder to detect in 
HI than their counterparts in low density regions. 

We thus conclude that, while the environmental depen- 
dence of the HIMF is only marginally statistically signif- 
icant, the evidence argues against morphological segre- 
gation as an adequate explanation of this dependence. 
It may be that the environmental dependence is caused 
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by increased deficiency in higher density regions towards 
smaller linear optical diameters. If this is the case, then 
it argues in favor of gas removal mechanisms that pref- 
erentially affect smaller galaxies. 

We caution, however, that this conclusion depends on 
the premise that galaxies of a particular morphological 
type would have the same HIMF shape within high and 
low density regions were we to exclude the effects of defi- 
ciency Croton et al. (2004) shows that the optical lumi- 
nosity function of the 2dF Galaxy Redshift Survey varies 
with environment even when one excludes the effects of 
morphological segregation. There may be similar varia- 
tion within the morphological subsamplcs of our data set 
that occur independently of the gas stripping that one as- 
sociates with deficiency. The statistics of our sample are 
insufficient to carry this line of inquiry further by segre- 
gating by morphological type more finely or segregating 
by morphological type and environment simultaneously. 

6. CONCLUSIONS 

We have extracted an optical diameter and HI flux- 
limitcd sample of 2771 galaxies from our HI archive to 
measure the HIMF of the local universe. The best-fit 
Schechter parameters for this sample are a = —1.24, 
log(M»/M ) = 9.99, 0, = 3.2 x lO" 3 Mpc" 3 . While 
the values of those parameters are as consistent with the 
Schechter parameters derived from HI blind survey sam- 
ples as the parameters from these samples are with one 
another, our estimate of the number of low mass objects 
is lower than the number predicted by most of these other 
surveys by a factor of a few. Because the low mass ob- 
jects are all in the very nearby universe where distance 
uncertainties are greatest and environmental differences 
in different regions of the sky are likely to be most pro- 
nounced, it is unclear to what extent this discrepancy 
is caused by environmental effects, systematic errors in 
distance estimates, poor sampling of galaxies with small 
linear diameters, or the exclusion of galaxies that should 
meet our optical diameter limit, but are missed by optical 
catalogs because of low surface brightness. 

Additionally, we have derived HIMFs for four mor- 
phological subsets of this sample. The HIMF appears 
to have a strong morphological dependence, with earlier 
types having flatter low-mass ends and higher values of 
M* than later types, a result which matches that of Z03. 
Our results for types Sa-Sc also appear to be roughly con- 
sistent with SGH, who found a declining low mass slope 
for those morphological types, despite that sample's poor 
sampling of the low mass end of the HIMF. 

We have also performed the first comparison of the en- 
vironmental dependence of the HIMF that compares all 
higher density regions in the survey volume to all lower 
density regions rather than restricting the comparison to 
a particular cluster or void vs. the field. Our subsamplcs 
have been divided by local matter density as determined 
by the PSCz density reconstruction. The two subsam- 
ples comprised of galaxies in higher density regions have 
HIMFs with flatter low-mass ends and lower values of 



M* than the subsamplc of galaxies in the lowest-density 
environments. We have argued that the flattening of the 
HIMF in the two higher density regimes is unlikely to be 
caused by morphological segregation, and that if it is a 
result of increased HI deficiency relative to the low den- 
sity regime, then this would be an indication of increas- 
ing deficiency towards smaller linear optical diameters at 
the low mass end of the HIMF, a result consistent with 
HI deficiency in clusters being caused by galaxy-ICM in- 
teractions. We cannot, however, rule out the possibility 
that the environmental dependence is caused by varia- 
tion within the morphological subsamples independent 
of the gas stripping in high density regions commonly 
associated with deficiency. The caveat involved in any 
conclusions that we may draw from the differences in 
HIMF shape among these three subsamples is that the 
differences are not particularly large, and in fact the three 
HIMFs are marginally consistent at the la level. 

The best way to resolve the question of how the HIMF 
depends on environment is with better statistics. If we 
had a sample with many more galaxies, we could shrink 
the error contours in Figure 10, and we could determine 
which galaxies are responsible for the variation in HIMF 
shape by segregating by morphological type and environ- 
ment at the same time, just as Croton et al. (2004) does 
with the optical luminosity function for the 2dF Survey. 
In this regard, upcoming HI surveys with the Arecibo 
L-band Feed Array may be able to resolve the issue. 

Better statistics may also be sufficient to resolve the 
discrepancy between our estimate of the low mass end 
of the HIMF and the values derived from most blind HI 
surveys. With a sufficiently large Hl-selectcd sample, one 
could compute the HIMF both including and excluding 
galaxies that were not a priori optically identified. It 
could then be determined whether our estimate of the 
low mass end of the HIMF is depressed by the exclusion 
of low surface brightness galaxies not included in optical 
catalogs. Likewise, better statistics would help in that 
they would improve our understanding of the environ- 
mental dependence, which could allow us to determine 
whether the low mass end of the HIMF is affected by 
local environmental effects. If the discrepancy in the low 
mass end of the HIMF is a result of errors in the flow 
model, however, then a larger sample is insufficient to 
resolve the issue. We would require either a better flow 
model or we would need a sample of low mass galaxies 
that goes deeper in redshift. 
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